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We consider the effect of inhomogeneous neutrino degeneracy on Big Bang nucleosynthesis for the 
case where the distribution of neutrino chemical potentials is given by a Gaussian. The chemical 
potential fluctuations are taken to be isocurvature, so that only inhomogeneities in the electron 
chemical potential are relevant. Then the final element abundances are a function only of the 
baryon-photon ratio r\, the effective number of additional neutrinos AN V , the mean electron neutrino 
degeneracy parameter £, and the rms fluctuation of the degeneracy parameter, cr^. We find that for 
fixed r), A7V„, and f, the abundances of 4 He, D, and 7 Li are, in general, increasing functions of ctj. 
Hence, the effect of adding a Gaussian distribution for the electron neutrino degeneracy parameter 
is to decrease the allowed range for rj. We show that this result can be generalized to a wide variety 
of distributions for £. 
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I. INTRODUCTION 

Many modifications to the standard model of Big Bang 
nucleosynthesis (BBN) have been explored (I] . One of the 
most exhaustively investigated variations on the standard 
model is neutrino degeneracy, in which each type of neu- 
trino is allowed to have a non-zero chemical potential ^ , 
and a number of models have been proposed to produce 
a large lepton degeneracy [||- ||. More recently, ob- 
servations of the cosmic microwave background (CMB) 
fluctuations have been combined with BBN to further 
constrain the neutrino chemical potentials ||- pl| . 

An interesting variation on these models is the pos- 
sibilitythat the neutrino degeneracy is inhomogeneous 
|2[- (l| . The consequences of inhomogeous neutrino 
degeneracy for BBN were examined by Dolgov and Pagel 
p3[ and Whitmire and Scherrer |jl5). Dolgov and Pagel 
examined models in which the length scale of the inhomo- 
geneity was sufficiently large to produce an inhomogene- 
ity in the presently-observed abundances of the elements 
produced in BBN. Whitmire and Scherrer investigated 
inhomogeneities in the neutrino degeneracy on smaller 
scales; in these models the element abundances mix to 
produce a homogeneous final element distribution. Us- 
ing a linear programming technique, they derived upper 
and lower bounds on the baryon-to-photon ratio rj for ar- 
bitrary distributions of the neutrino chemical potentials 
and showed that the upper bound on 77 could be con- 
siderably relaxed. However, the resulting distributions 
for the neutrino chemical potentials were quite unnatu- 
ral. Hence, in this paper, we examine a more restricted 
class of models, in which the distribution of the chemical 
potentials is taken to be a Gaussian. 

In the next section, we discuss our model for inhomo- 
geneous neutrino degeneracy. We calculate the effect of 
these inhomogeneities on the final element abundances 
and discuss our results in Sec. 3. We find that, in most 
cases, the effect of Gaussian inhomogeneities in the elec- 
tron neutrino chemical potential is to increase the abun- 
dances of deuterium, 4 He, and 7 Li relative to their abun- 
dances in models with homogeneous neutrino degeneracy. 



II. MODEL FOR INHOMOGENEOUS NEUTRINO 
DEGENERACY 

We first consider the case of homogeneous neutrino 
degeneracy. For this case, each type of neutrino is char- 
acterized by a chemical potential fii (i = e, /i, r), which 
redshifts as the temperature, so it is useful to define the 
constant quantity £j = i-ii/Ti. Then the neutrino and 
antineutrino number densities are functions of 
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and the total energy density of the neutrinos and an- 
tineutrinos is 
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Electron neutrino degeneracy changes the n <-> p weak 
rates through the number densities given in equations 
(|l|) and (g), while the change in the expansion rate due 
to the altered energy density in equation (J3J) affects BBN 
for degeneracy of any of the three types of neutrinos. (See 
Ref. M for a more detailed discussion). 

Now consider the effect of inhomogeneities in the neu- 
trino chemical potential. As noted in Ref. p5| , neutrino 
free-streaming will erase any fluctuations on length scales 
smaller than the horizon at any given time. Thus, in or- 
der for inhomogeneities to affect BBN, they must be non- 
negligible on scales larger than the horizon scale at n «-> p 
freeze-out, which corresponds to a comoving scale ~ 100 
pc today. On the other hand, if the neutrino chemical 
potential is inhomogeneous on scales larger than the ele- 
ment diffusion scale, estimated in Ref. to correspond 
to a comoving length ~ 1 Mpc, then the result will be an 
inhomogeneous distribution of observed element abun- 
dances today (the possibility considered in Ref. p3[ ). 

To make any further progress, we need a specific distri- 
bution for neutrino chemical potentials. In analogy 
with the distribution of primordial density perturbations 
(and in accordance with the central limit theorem) we 
take this distribution to be a multivariate Gaussian. Such 
a distribution is entirely characterized by the power spec- 
trum of fluctuations, P(k). For a power spectrum of the 
form P(k) oc k n the rms fluctuation on a given length 
scale A is given by 

a e cx A-("+ 3 ^ 2 (4) 

We wish to consider only cases for which the presently- 
observed element distribution (determined by at a co- 
moving scale of ~ 1 Mpc) is homogenous, while the dis- 
tribution is highly inhomogenous on the horizon scale at 
nucleosynthesis (a comoving scale of ~ 100 pc). Since 
our two length scales of interest differ by a factor of 10 4 , 
this condition can be satisfied for any k n power spectrum 
with n > —3. For instance, for a white-noise power spec- 
trum, n — 0, a value of er^ = 1 at the BBN horizon scale 
corresponds to er^ = 10~ 6 at the element diffusion scale. 

Given these conditions, it is a good approximation to 
assume that BBN takes place in separate horizon vol- 
umes, with the value of £ taken to be homogeneous within 
each volume. At late times, the elements produced within 
each volume mix uniformly to produce the observed ele- 
ment abundances today. 

We make the additional assumption that the neutrino 
fluctuations are isocurvature, so that the total fluctua- 
tion in energy density is zero, even when the chemical 
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potential is inhomogeneous. This implies that the over- 
density in the degenerate neutrinos is compensated by an 
underdensity in some other component. In Ref. O], for 
example, the degeneracies in each of the three neutrinos 
are arranged so that the total density remains uniform. 
In Ref. ]l4l i the compensation is produced by a ster- 
ile neutrino. Such models have the advantage that they 
produce no additional inhomogeneities in the cosmic mi- 
crowave background as long as the compensating energy 
density does not include photons or baryons. (Note that 
this is not the assumption made in Ref. [pH]). We also 
assume for simplicity that rj remains uniform in the pres- 
ence of an inhomogeneous lepton distribution. With this 
set of assumptions, the only neutrino for which inhomo- 
geneities in the chemical potential are important for BBN 
is the electron neutrino; the effect of the other neutrino 
chemical potentials is to alter the total energy density, 
which is now assumed to be homogeneous. 

It has recently been noted that if the large mixing 
angle solution of the solar neutrino problem is correct, 
then neutrino flavor oscillations will cause the neutrino 
chemical potentials to equilibrate prior to Big Bang nu- 
cleosynthesis []T6|-|l8||. In our inhomogeneous model, the 
effect of this equilibration would depend on the compen- 
sation mechanism for the inhomogeneities. In models in 
which the fluctuations in the electron neutrino chemical 
potential are compensated by fluctuations in the chemical 
potentials of the fj, and r neutrinos, the effect of such fla- 
vor oscillations would be to erase any spatial fluctuations 
in the chemical potentials. In models where the electron 
neutrino chemical potential fluctuations are compensated 
in some other way, the chemical potentials of all three 
species would be equal at any point in space, but the 
spatial fluctuations would be preserved. Any large AiVj, 
in this case would have to be due to some other form of 
energy beyond the standard three neutrinos. 



III. CALCULATIONS AND DISCUSSION 

The model described in the previous section can be 
completely specified by two parameters, the (inhomoge- 
neous) electron neutrino degeneracy parameter, £ e , and 
the additional (homogeneous) energy density due to the 
degeneracy of all three neutrinos plus any additional rcla- 
tivistic component. We parametrize the latter in terms of 
AN U , the effective number of additional neutrinos. This 
second parameter hides our ignorance about the compen- 
sation mechanism and about the degeneracies among the 
other two types of neutrinos. In our simulation, we take 
£ e to be homogeneous within a given horizon volume dur- 
ing nucleosynthesis. Different horizon volumes may have 
different values of £ e , which are given by the distribution 
function /(£ e ), i.e, the probability that a given horizon 
volume has a value of £ e between £ e and £ e + dt; e . (Since 
we are considering only inhomogeneities in electron neu- 
trinos, we now drop the e subscript). We take /(£) to 



have a Gaussian distribution with mean £ and rms fluc- 
tuation ctj: 

/(0 = -=^exp[-(C-a 2 /2a|] (5) 

Then the final primordial clement abundances, for a fixed 
value of 77 and AN V , will be functions of £ and tr^; we can 
write, for a given nuclide A, 

Xa = ^ X A (£)f(t)d£, (6) 

J — 00 

where Xa(0 is the mass fraction of A as a function of 
£, and Xa is the mass fraction of A averaged over all 
space; after the matter is thoroughly mixed, Xa will be 
the final primordial mass fraction. 

A full treatment for all possible values of 77, AN„, £, 
and 0£ is impractical. We have chosen to concentrate 
on variations in the latter two quantities, since we are 
most interested in the effects of inhomogeneities in the 
chemical potential. Because of the large number of free 
parameters and the difficulty of exhaustively searching 
all of parameter space, our goal is to discern any general 
results which are independent of r\ and AN V . 

There are now strong limits on r\ from the cosmic mi- 
crowave background alone, independent of BBN. We ex- 
amine two extreme values for 77: r\ — 4 x 10~ 10 and 
rj = 1 x 10 -9 ; these represent very conservative lower 
and upper bounds on rj from the CMB in models with 
non-zero neutrino degeneracy |icf| . For AN V , we con- 
sider AJV„ = and 5. Note that the first of these is 
only possible if the extra energy density in the degener- 
ate electron neutrinos is compensated by a decrease in 
the energy density in some other relativistic component. 
For each of these cases, we calculate the abundances of 
4 He, D, and 7 Li as a function of for £ = — 1 to +1 
in steps of 0.5. Our results are displayed in Figures 1-3. 
In each of these figures, we also show observational lim- 
its on the primordial element abundances from Ref. [fl9| : 
2.9 x 1CT 5 < (D/H) < 4.0 x I0~ 5 , 0.228 < Y P < 0.248, 
and -9.9 < log( 7 Li/H) < -9.7. 

The general behavior of the element abundances in 
Figs. 1-3 is very clear. As expected, for er^ <C £, 
the abundances of deuterium, 4 Hc, and 7 Li are un- 
changed from their values in the corresponding homo- 
geneous model with the same value of £ . At the opposite 
limit, when ^> £, the models all converge to a single 
limiting value; again, this is what one would naively ex- 
pect. What is interesting is that, with a few exceptions, 
the introduction of a Gaussian distribution of values for 
£ results in an increase in the abundance of each element 
relative to the corresponding homogeneous model with 
the same value of £. The only exceptions occur for 4 He 
with negative values of £, (for which Yp is far too large 
to be physically reasonable), and some of the 7 Li curves, 
for which there is a tiny decrease in the 7 Li abundance 
over a short range of values. 
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This result may seem surprising, but it is a simple con- 
sequence of the behavior of Xa(Q- In particular, if Xa(£) 
is a convex function (Xl(£) > 0), then Jensen's inequal- 
ity @ gives 

/OO 
*4£)/(£K > X A (0- (7) 
-oo 

We find, for example, for AiV^ = 0, and both values 
of ?7, that our X(£) curves are all convex in the range 
—2 < £ < 2, with the exception of 4 He at £ < —1, and 
7 Li with = 1 x 10 ~ 9 . These are precisely the regimes 
for which we observe equation (j^) to fail. Of course, 
none of the Xa (£) curves is convex for all values of £ ; the 
practical condition for equation (Q) to hold is that the 
X(t;) curves be convex as long as /(£) is non-negligible. 
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FIG. 1. The primordial 4 He mass fraction, Yp, as a function of the rms fluctuation in the electron chemical potential a^, for 
the indicated value of the mean electron neutrino chemical potential £. Each figure corresponds to the indicated value of the 
baryon- photon ratio r\ and the effective number of extra neutrinos AiV„. The gray shaded region gives observational limits on 
Y P from Rcf. [19]. 
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FIG. 2. As Fig. 1, for the primordial ratio of deuterium to hydrogen, (D/H). 
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FIG. 3. As Fig. 1, for the primordial ratio of 7 Li to hydrogen, ( 7 Li/H). 



This simple behavior allows us to draw some useful 
general conclusions. In models in which £ and AN V are 
allowed to vary freely, if we fix £ and trace out the allowed 
region in the r], AN U plane, then the upper and lower 
bounds on 77 are set primarily by the upper observational 
bound on 7 Li and the upper observational limit on D, 
respectively, with the 4 He limits serving primarily to set 
the bounds on AN V However, our results indicate 

that the general effect of going from a homogeneous to 
an inhomogeneous distribution in £ is to increase both 
the deuterium and the 7 Li abundances. (Again, we note 
a slight decrease in 7 Li over a small range in , but this 
is a tiny effect). Hence, the net effect of introducing this 
inhomogeneity will be to decrease the allowed range for 
r], in comparison with the corresponding homogeneous 
model. This is a rare example in the study of BBN in 
which the introduction of an extra degree of freedom does 
nothing to increase the allowed range for 77. Instead, the 
effect of adding a Gaussian distribution of values for £ is 
to decrease the allowed range for 77. 

Although we have assumed a Gaussian distribution for 
£, our results are much more general. In particular, as 
long as our distribution /(£) is negligible over the range 
of values of £ for which Xa(£) is not a convex function, 
we expect equation (^) to hold. This would apply, for 
example, to a top hat distribution with the same values 
of <T£ as those examined here. Moreover, the distribution 
/(£) ncc d not even be symmetric for our results to apply. 

Our results contrast with those of Ref. |^5J, which 
found an expanded upper limit on 77 in models with in- 
homogeneous £. The reason for this difference is that the 
models examined in Ref. |Te|| allowed for an arbitrary 
distribution in £, and large increases in rj occurred for 
bizarre distributions in £. In particular, the distributions 
in Ref. |D| sampled extreme values for £, outside the 
range for which all of the Xa(£) functions are convex. 
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